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Non-orthogonal multiple access (NOMA) with advantages such as superior 
spectral efficiency has been considered as a promising multiple access 
technique for the fifth-generation (5G) mobile networks. In this research, 
we propose energy harvesting (EH) NOMA system relaying network over 
Rayleigh fading channel using selection combining at the receiver. 
Eirstly, we investigate the system performance in terms of the closed-form 
expression of the outage probabihty (OP). Here we compare the OP of two 
destination users of the proposed system. Einally, all the results is convinced 
by the Monte Carlo simulation. Prom the results, we can confirm that all 
the analytical and simulation results are the same in connection with 
the primary system parameters. 
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1. INTRODUCTION 

Non-orthogonal multiple access (NOMA) with advantages such as superior spectral efficiency has 
been considered as a promising multiple access technique for the fifth-generation (5G) mobile 
networks [1-6]. In comparison with the traditional water-filling power allocation strategy, NOMA can 
transfermore power to the users in the worse channel conditions. It leads to a better trade-off between system 
throughput and user fairness. From the previous researches, authors investigated the impact of user pairing o n 
downlink NOMA systems as in [7], the authors in [8] have considered the power allocation with the max-min 
fairness criterion of the model system. An uplink NOMA scheme with joint power and subcarrier allocations 
has been proposed and investigated in [9]. In [10], a cooperation-based NOMA scheme for coordinated direct 
and relay transmissions is studied. Furthermore, a diversity-oriented detection mechanism for the cooperative 
relaying system using NOMA is considered by authors in [11]. The performance of transmit antenna 
selection for NOMA assisted multiple-input-multiple-output (MIMO) relay networks is presented in [12]. 
Then inspired by user collaboration, a cooperative NOMA transmission scheme is presented in [13]. 

In this research, we propose energy harvesting (EH) NOMA network system over the Rayleigh 
fading channel using selection combining at the receiver. In the first stage, we investigate and derive 
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the closed-form expression of the outage probability (OP) of the model system by Monte Carlo simulation. 
Then the comparison of the OP of two destination users is investigated in connection with the primary system 
parameters. EinaUy, aU the results are convinced by the Monte Carlo simulation. Prom the results, we can 
convince that all the analytical and simulation results are the same in connection with the primary system 
parameters. The main contributions of this research are listed as: 

The EH NOMA network system over the Rayleigh fading channel using selection combining 
at the receiver is proposed. 

- The closed-form expression of the system OP is investigated and derived. 

- The comparison of the OP of two destination users is investigated. 

AU the results are convinced by the Monte Carlo simulation. 


2. SYSTEM MODEL 

The EH NOMA network system over the Rayleigh fading channel using selection combining 
at the receiver is drawn in Pigure 1. In this system model, we denote Source is S, Relay is R, Di and D 2 
are two Destination nodes. We assume that aU links between them (i.e., S-to-Di, 2 , S-to-R, and R-to-Di, 2 ) 
are available. The Rayleigh fading channel coefficients of S-to-Di, 2 , S-to-R, and R-to-Di ,2 links are denoted 
as hsDi, 2 , hsR, and hRDi, 2 , respectively as shown in Pigure 1. The energy harvesting and information 
processing (IT) are iUustrated in Pigure 2. In this protocol, the transmission block time T, which consists 
of three-time slots. In the first-time slot aT (a is the time switching factor, 0<a<l), the relay harvests energy 
from the source node S. In the second interval time (l-a)T/2, the source S transfers the information to 
R and Di ,2 at the same time. EinaUy, the relay node R transfers the information to the destination nodes 


T 


aT (l-a)T/2 (l-a)T/2 


Pigure 2. The time switching protocol 


EH atR 


IT 

S^R, Dl, D2 


IT 

R^D1,D2 


Di ,2 in the remaining time slot (l-a)T/2. 



Pigure 1. System model 


Energy harvesting phase 

In the first slot time aT, we formulate the received signal as: 

yr 

The received energy can be computed as: 

E^=riaT\KfPs 

Where 0 < 77 < 1 is the energy conversion efficiency. 

The average transmit power at the R node can be calculated as: 


(l-a)r/2 


--K 


KfPs 


Where we denote k — 


Irja 

\-a 


( 1 ) 


( 2 ) 


( 3 ) 
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Information transmission 

In the second slot time, S transmits the superposed signal +^a 2 P^X 2 to relay, Di and D 2 

according to the NOMA scheme [1]. ai, a 2 are the power allocation coefficients for Di and D 2 , where we 
assume ay-\-a 2 =1,02 > ay,xi are the transmit signal for Di and D 2 , respectively. Hence, the received 

signal at the R node, Di andD 2 can be expressed as respectively [14]. 


= K + 4(hPs -*2 ) + 

(4) 

yi\ ~ Kd^ ^l^2^s -^2 ) ■*■ 

(5) 


(6) 


Where additive white Gaussian noise (AWGN) at theDi and D 2 with zero-mean and variance 

No, respectively. From (6), we separate the received signal at the D 2 into two parts included signal and noise 
part for D 2 to detect x 1 . 

= Kd, + K, ( 7 ) 

signal noise 


Hence, the received signal to interference and noise ratio (SINR) for D 2 to detect xi is given as: 



Isignal^ 

1_ I' 

h 

sd2 

f 


|2 

e| 

f| . |2 
^\noise\ 

r 




Q^ij/ +1 


( 8 ) 


Where y/ = is the transmit signal to noise ratio (SNR). 

^0 

From (5), SIC is first used for Di by detecting and decoding the information of D 2 . 
Hence, the received signal at the Di can be rewritten as: 

Ja = V ^/^^2 + K, \ (9) 

signal noise 


Then, the received SINR atDi is given by: 

ysd 2 i “ I ,2 

\Ka\ 


( 10 ) 


After the far user message is decoded, Di can decode its own information with the SINR 
as following 


Ysd, 



( 11 ) 


During the third slot time, the relaying node will amplify the received signal and forward to Di and 


D 2 with the amplifying factor (3 = — = I- ^ -. The received signals at the Di and D 2 

are expressed as,respectively. 
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yrd, = Kd, Xr + = Kd,Pyr + 


Substituting (4) into (12) 

yrd, = Kd, P K )+"'■]+"'■rfi 


And 


yrd, =Kd,Xr+nrd, =Kd,Pyr+n^d, 

= Kd, KPy^^’^ + ^/^•^ ) + Kd, P^r + 


( 12 ) 


(13) 


(14) 


Where P^rdy’>^rd^ was denoted the AWGN at Di and D 2 , respectively. Erom (14), the received SINR for D 2 
to detectX 2 is given as: 



\signalf 

!_ \h 

|2 
sr \ 

^rd) 

2 

a^y/ 

Ej 

M . |2) 

\noise\ j 

1 

1 . 

\¥+Kd, 

^ +\i p 


(15) 


Using the fact thatNo«Pr, hence, (15) can be rewritten as: 

Yrd — -^- - - 

\Kt\Kd] aiW 

\Kt\Kd,\ +1! K 

The received SINR when Di first detect D 2 ’s information is given by: 

\K'\\Kd!((h¥ 

\Kf\Kd!\(h¥+\Kd!\+^if^ 

Applying SIC operation,the received SINR for Di is given by: 

\Kt\Kd!\ (h¥ 

^rdi ^ |2 

Kd] +1/K- 


(16) 


(17) 


(18) 


3. SYSTEM PERFORMANCE ANALYSIS 

The outage probability (OP) of the destination node D 2 can be given by: 

= Pr [max , y ^,^) < ^^02 ] 

=^(rsd, <ro 2 )Pr{rrd, <^ 02 ) 

Where Yq 2 — “1 with R 2 being the target rate at the D 2 . Prom (8), we have: 
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^{Tsd, < ro2 j = Pr 7^12- < ^02 

a,i,r + \ j 

= Pr f ¥ [«2 - ro2«i ] < ro2) 

, { Yta^sd, 1 -x ^ 

1-exp-f-^T , ifa 2 >YfaCh 

y ¥[a 2 -M\) 

1 , ifa^<Y^^a^ 

Where ^Sd2 is the mean of the random variable (RV) . From (16), we can calculate the second 

probability of (19) 

Pr{r.rf, <^ 02 ) 

\Kt\Kd.ta^¥ 

- Pr -- - -2- < 7 w. - 

a,y/+\h^,\ +\Ik ^ 

J -M 2 E-= 

(pM 2 a.ii/ + (p 2 +\l K 

\ ^ 

Pr (p2{(Pia2i/^-Yo2¥ai¥-ro2)<— = 

K 


]fnd(p, + \^\(P2<—r -^^ 

i i y n¥ia2¥-ro2¥iai¥-ro2\) 

1 - 4 J exp I—f--.1 exp(-4^(?, )d(p, 

i y ’^ma2¥-ro2¥ai¥-ro2\) 


Where we denote ^2 ,(p^-\h ,(p 2 =\h^ Ar ^rd 2 are the mean of RV respectively. 

By changing a variable -ro 2 ^,a,i^-ro 2 ^<p, = > (21) canbe rewritten as: 

r 02 


^ ^ ro 2 

f ( ^024,) ( 

exp --^ exp — ax 

0 y ) I /o 2 ) 


Apply eq (3.324,1) of [15], (22) canbe reformulated as: 


'■{/rd, <ro 2 ) = l-2exp(-45) 


ArAc/ ^ ^sr^rd ^ 

^ rd2 ^ sr_j^ 


Where K^{^) is the modified Bessel function of the second kind and v^^ order. Substituting (20) and (23) into 
(19), finally, we can obtain: 


OPd^ ={l-exp(-4'9)} 


-2exp(-4.9)x. 


/44 5 


-^>cKA 2A-^ 
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Hence, the outage probability (OP) of the destination Di can be given by: 
<ro2^nrf, <roi) 

<ro2’rrd, <rm) 

=[i-Pr{nrf„ ^roi)]x 

[l-Pr (/rd,, ^ rm . Trd, ^rm)] = Pi ^Pi 


Where =[l-Pr(r,,^, >ro 2 .n., ^Xoi)]’ ^ro.)]' ^oi =2"'"'-1 with Ri being the target 

rate at the Di Erom (10) and (11), Pi can be calculated as: 


Pi- 


l-Pr 




hsd, 


^ria’\hsd] m^rrn 


- 1 -Pr 




Via 




' ' a^y/) 


= Pr 


max 


r<yi 

Ym 11 

_api/-r^a,y/' 

'a,y/\\ 


|4rf,f<S =l-exp{-/l,^s) 


(26) 


Where is the mean of RV |/z^ | and 
calculated as: 


^ = inax 


ro2 roi 


. Prom (17) and (18), P 2 can be 


P2=l-Pr 


\K\^ Kd, 


\Kf\Kd,\ m+\Kd,\ +1/; 


- 7o2 


\hJ\KdSm. 

Kd^ + 1 /'^' 


Toi 


= 1-Pr<i(33 >A 

(Pi-A 


R I K- 

^P^\(p3 ^- -,<Pi 

(Pi-B 


= 1-Pr<j^(?3 >2=^,^ >S 

(Pi-- 


00 00 


= >-( 
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exp 
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exr 


(27) 
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(p^=\Kt^%=\Kd\ > 


Where 


^02 


,B = -^,S = max(A,B) 


2 I 1^ 

.And are the mean of RV \h^^\ , 


respectively changing variable y = ^—S, (27) can be rewritten as: 


^ 2 =l- 2 j 


exp - 

0 I 


lexp[-4(y-HS)}fy 


°° f A/ 

= 1-4 exp(-4^)1 exp 

0 I 'e}- y 


(28) 


dy 


Apply equation (3.324,1) of [table of integral], (22) can be reformulated as: 


Pj =l-2exp(-2,^5)> 


42re(|“ 


K 


xKi 


2 l^srAdi^ 


(29) 


where K^(^) is the modified Bessel function of the second kind and v^^ order. Substituting (26) and (29) 
into (25), finally, the outage probability ofDi can be claimed as: 


=[l-exp(-4H)] 
l-2exp(-4H)x 






U 

2 I sr'^rdi 


K 


(30) 


4. NUMERICAL RESULTS AND DISCUSSION 

In this section, we investigate the OP of the model system using Monte Carlo simulation in 
connection with the primary system parameters [16-26]. In Figure 3, the effect of a 2 on the system OP 
is plotted with the primary system parameters as r|=0.8, Ri=R2=0.25, 0.5 bps/Hz, a=0.5. From the results, 
we can see that the OP of the destination Di significantly decreases with the rising of a 2 . However, the OP 
of the D 2 decreases when a 2 increases from 0.55 to 0.7 and after that has a massive increase with 
the remaining values of a 2 . In the same way, the system OP versus a is drawn in Figure 4. In this simulation, 
we set ri=0.8, Ri=R2=0.5 bps/Hz, Ps/No=0.5, and a2=0.7, 0.9. As shown in Figure 4, we can state that the OP 
of the model system has a slight decrease when a varies from 0 to 1. In both Figures 3 and 4, the simulation 
andanalyticalresults are the same with all values of a anda 2 . 
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Furthermore, we investigate the effect of rj and \|/ on the OP of the model system as plotted in 
Eigures 5 and 6. In these Eigures, we set the primary parameters as a2=0.7, 0.9, Ri=R2=0.5 bps/Hz, 
Ps/No=0.5. Prom Eigures 5 and 6, it can be observed that the system OP has a slight decrease when p varies 
from 0 to 1 and the system OP crucially decreases with the rising of \|/ from 0 to 25. In all two Figures, 
the analytical and simulation results agree well with each other. On another hand, the comparison 
of the system OP of two destination nodes is demonstrated in all Figures. From the results, we can state that 
the system OP of the destination node D 2 is better than the destination node Di. 




Figure 5. OP versus p 


Figure 6. OP versus \|/ 


5. CONCLUSION 

In this research, we propose EH NOMA network system over the Rayleigh fading channel using 
selection combining at the receiver. The closed-form expression of the OP of the model system by Monte 
Carlo simulation is derived. Moreover, the comparison of the OP of two destination users is investigated in 
connection with the primary system parameters. From the results, we can see that aU the analytical 
and simulation results are the same with the primary system parameters using Monte Carlo simulation. 
These results can be provided a novel recommendation for improving the performance of the EH NOMA 
communication network system. 
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